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Extended f{R) theories of gravity have been investigated from the symmetry point of 
view. We briefiy has been investigated Noether symmetry of two types of extended f(R) 
theories: f{R, T) theory, in which curvature is coupled non minimally to the trace of 
energy momentum tensor and mimetic f(R) gravity, a theory with a scalar field 
degree of freedom, but ghost-free and with internal conformal symmetry. In both cases 
we write point -like Lagrangian for fiat Friedmann-Lemaitre-Robertson-Walker (FLRW) 
cosmological background in the presence of ordinary matter. We have been shown that 
some classes of models existed with Noether symmetry in these viable extensions of f{R) 
gravity. As a motivated idea, we have been investigating the stability of the solutions and 
the bouncing and ACDM models using the Noether symmetries. We have been shown 
that in mimetic f{R) gravity bouncing and ACDM solutions are possible. Also a class 
of solutions with future singularities has been investigated. 

Keywords: Modified gravity theories; Cosmology; Noether symmetry; dynamical sys¬ 
tems. 


1. Introduction 

A challenge to the contemporary relativistic cosmology is provided by a set of obser¬ 
vational data, indicated on late time acceleration expansion of the whole Universe 
as well as the initial era, the inflationary epoch ra-i- In the framework of the gen¬ 
eral relativity (GR), inflationary model, quasi stable de Sitter cannot be realized 
as a possible physically acceptable model without any extra scalar matter field. It 
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is needed to include some matter fields, like scalar field(s) to resolve it. Another 
approach which is totally revolutionary is to replace Einstein gravity by some ex¬ 
tended forms of classical geometrical objects. A reasonable candidate to explain 
this situation is modified gravity, in which we replace Einstein-Hilbert action, given 

by 

Seh = [ + SehIom ( 1 ) 

Jm 2k 

(Here SehIom is the boundary term) by another set of geometrical objects like sec¬ 
ond order invariants R^''°‘^ Ri^vap , ■■ or functions of Ricci scalar R, Gauss- 

Bonnet topological invariant G = R^''°‘d -1- R? and so on (see [S]- 

m for reviews). Historically, to be more precise, the simplest potentially reason¬ 
able candidate was f{R) gravity, a theory which it was proposed originally before 
recent activities m) and later motivated in light of the recent observational data, 
as a valid, physically reasonable, ghost-free and stable alternative theory instead of 

GR 0 - [g. 

Another type of modified gravity is the one in which geometry has been coupled 
to the matter fields non- minimally (see |16) for An updated review of such models). 
Different types of non-minimally coupled models have been proposed like f{R,Crn) 
where Cm stands for matter Lagrangian and f{R,T) [TT] .where T is the trace of 
the energy momentum tensor of matter fields,which is defined by 


T — — 


2 d —gLm) 




Or after a simple checking, it can be rewritten as the following: 

dLm 


Cfxi/ — g^i^Cm 2 




( 2 ) 

( 3 ) 


In this theory, T = T^. Because of simplicity and beauty form, this theory attracted 
several activities in literature |18j-[33). Symmetry is an important issue to be ad¬ 
dressed in any physical system under study. There are two classes of symmetries: 
global symmetries , in which the physical system( dynamical system) respects some 
types of transformations, which are defined by functions of coordinates. Another 
is local, in which the conservation law gives us a ’’local” conserved quantity like 
charge. 

As we mentioned it before, f{R) gravity is a ghost free, and conformally equiva¬ 
lent to the scalar field theory in Einstein frame. There are several interesting features 
in this theory to be useful for late time acceleration, dark energy and dark matter 
halo problems. Also, inflation can be realized successfully using this simple and ef¬ 
fective modified theory of gravity. Recently to resolve dark matter problem, a new 
type of modified gravities has been proposed as titled mimetic model and modified 


®‘This common definition of energy-momentum tensor can not be used to fix the form of the 
Lagrangian. For example, of perfect fiuid there is no way to read Lagrangian from the form of T^i, 
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versions of it [33]-[12]. Basic hidden idea behind this new modification of gravity is 
to propose a conformally invariant, scalar theory of gravity in which scalar degree 
of freedom does not cause any problem with ghosts. The way is to reparametrize 
metric as a conformal transformation of an auxiliary metric unphysical” metric. 
The point is, the scalar field appeared as conformal function and plays the role of 
an internal degree of freedom. This scalar field is not ghost and it is constructed 
to have unit norm = —1. This norm is defined on physical metric. It is 

remembering for us the role of velocity of a test particle with unit norm in the 
comoving frame of particle coordinates. An interesting feature is if we write FLRW 
cosmological equations, an extra term proportional to a“^, appeared. It mimics 
dark matter. So, it was adequate to name it as “mimetic dark matter” or briefly 
mimetic gravity. To unify f{R) gravity with this very interesting mimetic theory, 
it has been proposed mimetic f{R) gravity [43] as a new class of modified gravities 
with the same inspiration as mimetic theory. Because of its complexity and more 
physical solutions, this new mimetic f{R) deserves further physical investigations. 
Very recently the dynamical behavior of mimetic f (R) has been investigated [44] . 
Our aim here is to address Noether symmetry issue of such mimetic models. In 
literature Noether symmetry has been investigated for different types of modified 
gravity like /(R), /(T), Galileons and so on [46]-[60]. Our aim in this paper is to 
explore Noether symmetry for the mentioned above to modify gravities: mimetic 
/(R) and/(R,r). 

The present paper is organized as follows: formal framework for /(R, T) theory 
of gravity is presented in Section [3] The formalism of f (R) mimetic theory is 
motivated in Section [3] The Noether symmetry approach is well understood briefly 
in Section 13] Noether symmetry is applied to /(R, T) in Section [S] and for /(R) 
mimetic model in Section [3] We summarize in the last section. 

We adopted a ”God” given system of units G = c = 1, where the gravitational 
coupling constant is given by = Stt. 


2. Formal framework of f {R, T) gravity 


Let us to start by the following simple extension of /(R) gravity in four dimensional 
spacetime: 




(4) 


In the action, R stands for the Ricci scalar of the Riemannian space-time, T is 
trace of the energy, momentum tensor of matter Lagrangian Rm, is defined by the 
simple expression T = Trace is an important quantity in quantum theory 

of inflation and quantum gravity. The dynamical quantity is just metric So, 
there is only a single field equation for it. If we calculate the variation of action we 
obtain: 


6S = 


IGtt 


fR{R,T)6R- 


h («> r) («. T) sr + 


^dX&] 
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In this variation, we define : fii{R,T) = df{R,T)/dR and friR^T) = 
df {R, T) /dT, respectively. Since R is the only geometrical quantity in our the¬ 
ory, we recall standard variational expressions for it as the following: 

5R = 5 {g^'^R^.) = + g^'^ (Va<5F^, - V.JF^a) , (6) 

A simple check point is to evaluate it in a local flat coordinate frame, when F^^ = 0. 
A variation of the connection termination, F gives us: 

= \g^^ {V^Sg^c. + yjgc^. - , (7) 

Consequently, we can write the final simplified form of variation of R as the follow¬ 
ing: 


5R = R^Jg^’^ + g^.aSg^’^ - . (8) 

By plugging these expressions in the total variation of action, 5S, we find: 

[fR {R, T) R^.Sgi^'' + Jh {R, T) - Jr {R, T) 

+/t (K, t) (ii, T) V" + 16^^ 


^d\si) 


To simplify more these functional, we need to redefine an auxiliary tensor field 
as a part of the variation of trace 5T as the following: 

<5 {g^^^T^p) 


Sgf^'' 


= T, 




0 


fll' 5 


( 10 ) 


And in a similar form as we define T, we are able to define trace of as the 
following: 


0 = 

^ Sgf^- ■ 


( 11 ) 


Using these simplifications, we write the following form of the equation of motion 
of/(R,T): 

Jr {R, T) R^.-\f (R, T) g^, + {g^^U - V^V,) /r (R, T) = Stt^ - It [R, T) - h (R, T) 0^2) 

It is remarkable that, when geometry is decoupled from the matter part, in the 
limit f{R,T) = f{R), we recover equation of motion of f{R) gravity from (fT^ . 

Trace of (IT^ gives us the following equation: 


fR (A, T)R + SD/fl (i?, T) - 2/ [R, T) = 8ttT - fr [R, T)T - fr [R, T) 0 (13) 
An alternative form of m is obtained as the following: 

fn {R^T) - ^Rgij,^^ + lf{R,T)gf,„ = Stt - fr {R,T) 


-/T(R,r)(0^,--05^, 



+ V;,V,/«(i?,T) . 


(14) 
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Conservation of the energy, momentum tensor in this type of modified gravity is 
checked by computing the covariant divergence of Eq. dm. Using the following 
geometrical identity we obtain: 




fn (R, T) [R, T) - V^V,) /« (i?, T) 


= 0, (15) 


So, conservation of energy, momentum tensor, gives us the following vector current 
term: 

^ [{Tgo + Qyu) In fr (R, T) + . (16) 

For a simplicity, we consider the flat FLRW metric in the following form: 

ds^ = —dt^ + a’^{t) 5ijdx''dx\ (17) 


Here a{t) stands for the scale factor. If we write down Lagrangian of f{R,T) for 
this metric and if we assumed that the Universe is filled with matter fields with 
effective pressure p and energy density p, we obtain T = Zp — p. Since this equation 
is a constraint, and because of R = 6(| + (f)^), we introduce a pair of Lagrange 
multipliers as the following {A,p}. The point Like Lagrangian after an integration 
part by part is written as the following (we set 2k^ = 1): 

£(a, i?, T, d, R, T) = (/(i?, T) - RfaiR, T) - T/t(R, T)) (18) 

-6(ad^fRiR,T) + a^dRfRR{R,T) + a^dtfRT{R,T)) - p{a)a^ + fT^Ma) - p(a)) 

Where we suppose that Cm = —p{a)- If f{R,T) = f{R), the Lagrangian density is 
written: 


£(a,i?,d,i?) = a3(/(R) - i?/fl(R)) - 6(ad‘^ fR{R) + a^dRfRRiR)'^ - p{a)a%l9) 


This is the form of Lagrangian in f{R) theory. We suppose that (fTHl) is an accept¬ 
able assumption for many cosmological applications, like matter dominant era or 
radiation. The associated equations of motion are given by a set of second order 
ordinary differential equations, Euler-Lagrange (EL) equations, are given by the 
following: 


dt dq 



( 20 ) 


for q= {a, i?, T} they are obtained as the foloowing: 


^RIrr = —3/ -I- -I- 3T/t — 6 —Jr — l2—fRR — 12 —/rt + ap'{a) -I- 3p(a) — 9 /tp(^)) 

a a 

+3/rP (a) — ^afrp'ia) + afTp'{a) (a) — GR^fRRR — 12RTfRRT — ^T'^/rtt — ^/rtT — 

^fRR~ — fRR a ~ (3p(a) - p(a))'j (22) 

• 2 

(R - 6^ - 6-)/flT + /tt(T - 3p(a) + p(a)) = 0. 
a 


(23) 
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To pass the case of f{R,T) = f{R), we obtain: 

QRfRR = -3/ + SR/r - 6 ^fn - U—Jrr + ap'{a) + 3p (a) (24) 

a 

-^R^Irrr - 

a 

R = (25) 

\ a aJ 

Our aim in this paper is to investigate Noether symmetry issue of (HU. Briefly, 
we are interested to know how Noether symmetry is able to ’’fix” mathematical 
forms of {f{R,T),p{a)}. 

3. formalism of mimetic f{R) gravity 

Although f{R, T) provides a reasonable and good extension of f{R) theory, it does 
not respect conformal symmetry. Also, extra degrees of freedom are possible. So, 
instabilities due to ghosts probably are happening. To resolve conformal symmetry 
and to be ghost-free, a model recently proposed as titled Mimetic F{R) gravity [43]. 
It is inspired from the mimetic theory [34135136137139140142] , a model in which dark 
matter problem is resolved as an integration constant. Also, it is self consistent with 
conformal symmetry. The basis of any type of mimetic theory is to parameterize of 
the Riemannian metric tensor as the following conformally transformed formula 

m 

9/1^ = -g^'^dp(l)da(l)gf,u , (26) 

Here we introduced a pair of auxiliary objects: the first is an auxiliary metric 
(unphysical and without dynamics) g^i, and the second is a scalar field degree 
of freedom cj) which has generally ghost, freedom. It is well known that using an 
orthogonality of metric, this scalar field satisfies the following constraint equation 
of motion: 

gi.gp.vAY''= (27) 

If we know the background metric g {g^v, ttiis equation fixes the form of the 
scalar field. If we interpret d^4' tti® components of a four velocity u^, then nor¬ 
malized = — 1 implies a possible normalization of the 4>. This normalization 

can be understood as the first integral of the equation of motion for cj) . We em¬ 
phasize here that the auxiliary metric g^i, is an internal object of the space-time 
manifold Following [43] . we write the following action for mimetic f{R) gravity in 
metric formalism : 



Here = Stt, R is the Ricci scalar which is computed by the physical metric 
and we also include the matter Lagrangian by Cm ■ What we need is to parametrize 
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the physical metric according to (Ei) . We rewrite the action of theory in the 
following equivalent form: 




2^2 


Cn 


(29) 


Thus, we perform variation with respect to auxiliary metric , we obtain: 

{R {gtj.u,4>)) - R {gtiu, fn {R {g^j.^, </>)) 

+v [g V [g </')) 

-9 {9/1^: □ {9tJ.i^, 4>) Ir {R {gi^uA)) + 

+d^(t)d^4>\2f {R ig^,u,(l))) - R{g^y,(t>) fR{R{g^„,,4>)) 


-3n {g Ir (R '/>)) + 


= 0 , 


(30) 


As a convention, here fn means df{R)/dR, and □ are different derivative 
operators with respect to Also by using a similar ’’dictionary” as we used in 
/(i?,T), we define T^jy as the effective matter, energy-momentum tensor for Cm is 
given by o . Variation with respect to the scalar field (f) gives us : 


{g ig/j.i', 2/ {R {g^y, (j))) - R {g^^, cj)) Jr {R (g^^, (j))) 

-^o{9i9fj^i^,(l})^J) Ir {Rig (/})) + }=0, (3i) 


Here like f{R,T), we define trace of the energy-momentum tensor as T = 
g {gfj,!^■> ■ It has been proven that this new theory is conformally invari¬ 

ant and ghost free [43]. So, it is remarkable to consider it as a valid extension of 
f{R) gravities. 

We consider the same FLRW metric as dm. In this case, the constraint equation 
leads to (j) = t where t is cosmic time. It is not so hard task to write FLRW equations 
for motion. Actually, because our aim is to investigate Noether symmetry, so what 
we need is just point like Lagrangian of the mimetic f{R) scenario. To be more 
generally speaking, we slightly modify the original mimetic f{R) by including a 
potential term and by introducing a Lagrange multiplier A as the following: 


5 = 




"\J-gigtivA) 


f{R (ff/.*!/,'/’)) 


2k2 


- V((/)) - + !)+£„ 


• (32) 


The point like Lagrangian for a fluid with pressure p(a, d, (p) and by taking in to 
account that Cm = p{ol, d, <f), is written in the following form: 

C{a, R, ct>, d, R, = a^fiR) - a^RfR{R) - 6(adVfl(Al) + a^aRfRR{R)) (33) 
—p{a, d, 4>)a^ — V— A(1 — 4P')a?' 
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Note that in general A = X{t). Equations of motion are written in the following 
forms: 

= -QR^Irrr - adpaa - aapaa - apat - 3apd + aPa (34) 

~IrI ~ 3R + 6(—+ 12— — 3(/ — A + A^^ — V{(p) — p) \ — 12-— 

\ a a J a 

- = i(i?-6(-)2) (35) 

a o a 

^ + 0(y + -) + J^ = O. (36) 

A a IX 

The second equation is just the standard definition of Ricci scalar for ELRW metric. 
The third one is reduced to the Klein-Gordon equation in the case of A = X{t) = 
Constant. The first equation becomes familiar as the equation of motion in f{R) 
gravity, if we set p = p{a), A = 1, = t, = 0. Clearly this equation posses 

de Sitter solution as R = Rq, a{t) = aoe^°*. 

Because of the importance of this model, we will study fixed points of the 
associated dynamical system, corresponding to this last case , when p = p{a), p = t. 
We applied it before in the context of general relativity [IS]. Here we review the 
basic concepts of an non-autonomous system. 

Consider the following differential equation for an non-autonomous dynamical 
state vector x: 


x = f{t,x) (37) 

f : [0, oo) X n ^ (38) 

iA = {feK” |||f ||2<0} (39) 

The equilibrium point, or fixed point is located at cc = 0 if and only if it solves the 
following algebraic equation for an instant of time,namely t: 

/(t,0) = 0,vt>0 (40) 

We define the Jacobian matrix ; J = [^] must be bounded function of t on a 
finite domain 7) and furthermore it satisfies smoothly the Lipschit’z lemma , as the 
following: 

II - f(t,y) ||< L II f- 2 / 11 , II f ||p= (^ I |P)p,l < p < oo (41) 

i 

There is an important theorem about the asymptotic stability of the system in the 
vicinity of the equilibrium point: 

Theorem I : It is possible to linearize the system of equations in the vicinity of the 
fixed point in the following form: 

Sx = A(t)dx 


(42) 
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(43) 


Theorem II: Suppose that a; = 0 be the hxed point of the system x = f{t,x) , and 
it satisfies the following auxiliary conditions: 

/: [0,oo) xD— 

= {a; G K" III X || 2 < r}, 

It is adequate to define the time dependent function A{t) = 1^=° ■ 

that the system has an exponential stable equilibrium point of the linearized 
system (HU) , then this point is the exponential stable equilibrium point of the 
nonlinear system dSzl) 

In this case the system of equations reduces to the following form: 

^IrrR = —QR^Jrrr + aPa (44) 


aR 


-fR{ - 3i? + 6(^)2 + 12^ - 3(/(i?) - A + A02 _ -p))- 12^ 


“ = i(R-6(“)2) 

a o a 


(45) 

(46) 


The first attempt is done by rewriting the system of equations in terms of a dimen¬ 
sionless ’’time” coordinate N = logo and a new set of dimensionless parameters 


as X 


A 


= r'X = log A}. Equations read as the following 
non-autonomous system (due to the potential term V{(f>) = V{N): 


— {h — 


h' = 2{--h) 
h 

r'=X, 


QHpfr 

frr 


(47) 

(48) 

(49) 


fr 

Um^frr 


6r 4 r 3 

lih ~ ’ ~ 2H^ 


{f-ViN)-piN)) 


X 


72H^frr 


C' = -3 - 


V{N) e 


-C 


(50) 


V{N, 


a~Cc 


= -3. 


(51) 


2Hp h 

Stationary (fixed) points are located at : 

f{r,) = V{N,)+p{N,)=0, X, = 0, r, = hl 

ZUp tie 

Because of the critical point is function of N, it’s moving when time is running. 
The Jacobian of the linearized system is given by the following: 

"Oi 02 0 0 

0 0 10 
J = 

Cl C2 C3 0 
_di 0 0 d4_ 


(52) 
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Where 


Qi = 

Ci = 
di = 


dfs 
dq^' 
dfi 
dq^ ’ 


<?* = {h,r}, 
q" = {h,r,X}, 

q" = 


And: 


h = \ - h), 


h = - 


(L 6-^0/rrr \ 

\h frr ' 


2X 


frr 

fr r 6 r 4 , r 




144ij4/„|. h'^^h^h 2iJ2/i2(/ 72H§f, 


OJrr 


/4 = -3 - 


ViN) e 


-C 


2Ho h 

The characteristic equation is given by: 

+ (—(i 4 — C 3 — ui) A^ + ((i 4 C 3 + d^cLi — C 2 + 0301 ) A^ 

+ (c^4C2 — (^40301 — Cl 02 + C 2 ai) A — ^4 ( — 0102 + C 2 ai) = 0. 
one eigenvalue is Ai = ^ 4 . Another eigenvalues read as the following: 

-f 12 C 2 — 4 C 3 O 1 -f 4 C 3 ^ A 4 Oi^ A 2 C 3 A A 2 oi A 


A 2 — 


, 1 C 2 

A3 — - TT X - 

12 X 


6A 


1 C 3 O 1 1 C 3 ^ 1 Oi^ 

-q-Q-f q C 3 A - Ol 

3 X 3 X 3 3 


1 


1 


3 X 

1 
3 


(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 


1 /TT iv^C 2 Aiv^C30i Aiv^C3^ Aiv^oi^ 

-T:t*v3x-A^-^-^- 

12 XX XX 


A 4 = -A 


3 - 


(62) 


In the above expressions : 

A = 12 A A <S oj ^ cjoj ^ A (—cs^ — 72 cg) aj + 36 C2C3 + 8 cs^ + 108 cj o^, ( 63 ) 

E"^ = 54 C2C3C1O2 — 108 C20iCia2 — 18 03^010102 — 18 0301^0102 — 12 C2^ A 81 ci^02^ — 3 02^03^ ( 64 ) 

A 24 C2^oi^ — 12 C20i‘* A 6 cs^oi^ — 3 03^01^ — 3 03^01^ A 12 010203^ A 12 040201^ — 24 02^0301 

— 6c2C3^0i^ — 6 C2C3^0i A 24 c 2 C 30 i^ 

X^ = 36 C2C3 — 72 0201 — 12 03^01 — 12 0304^ A 108 0402 A 8 03^ A 8 oi^ A 12 p, ( 65 ) 

= 54 02030402 — 108 C2O1C1O2 — 18 03^040402 — 18 0304^0402 A 81 04^02^ — 3 02^03^ A 24 02^04^ (66) 

— 12 C20i'^ A 6 03^04^ — 3 cs'^oi^ — 3 cs^oi"* — 12 02^ A 12 040203^ A 12 040204^ — 24 02^0304 — 6 0203^04^ 
-60203^04 A 24 020301^. 
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To have stable solution we must have ^4 = —3 + '^ 2 Ho^ ^ Another eigenvalues 

Ai, 2 < i < 4 must satisfy 5RAi < 0. It is remarkable to mention here that the stable 
manifold is defineb by the de Sitter space-time. So, the model for time intervals 
N > Nc has stable de Sitter solution. 


4. Noether symmetry approach to dynamical systems 


Symmetry is an important issue to be addressed in any physical theory. The basic 
property of a system with a defenite type of symmetry is the existence of an asso¬ 
ciated conserved quantity under this kind of symmetry. Let us consider a typical 
dynamical system, is defined by a set of configurations coordinates set qi. Gener¬ 
ally speaking, the dimension of the system is defined as the number of independent 
coordinates of the system. We assume that the dynamics of the system are defined 
by the point like Lagrangian is given by:L = 1 < i < N. For each 

coordinate, it is possible to define a ’’unique” first order conjugate momentum: 



(67) 


Euler-Lagrange equation of motion is given by the following set of N-ordinary second 
order diffrential equations: 


■ n 

”•- 97 =“' 


( 68 ) 


What we call it as Noether Symmetry Approach^^-^\ is the existence of a 
’’unique” vector field X on tangent space TQ = {qi,qi}): : 


X = T, 


N 


a\q) 


dq^ 


\q) 


dq^- 


(69) 


If we can find ’’generators” coefficients ai{qj), then we can strictly say that our 
dynamical system must satisfy the following geometrical constraint (is called as Lie 
derivative of the Lagrangian : 


LxC = 0 


(70) 


where explicitly we have: 


dC 


dC 


LxC = XC = El, a\q):^ + a\q)^.^^ 


or equivalently we can write it as the following: 


-^N 


dt 


def 


= LxC. 


(71) 


(72) 


It is an easy task to show that, existence of Noether symmetry implies that the 
system has the following conserved (local) quantity: 

V~' \^N i 

ZjQ — Pi 


(73) 
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The equations are obtained setting to zero the coefficients of the terms qi'^q‘j,0 < 
m + n < N in GB). They are several interesting applications of this symmetry 
approach in different cosmological models in different models |50)-|60|. In our paper 
we’ll apply this approach to f{R,T) and mimetic f{R) theories of gravity. 


5. Noether symmetry in f{R, T) gravity 

The plan in this section is to study the system of Noether equations for Lagrangian 
given by (fTOl) . If we write down (ITT]) equation for (ITOl) , 


dC OL / .Oa ■ tta ^ e'er \ cH- 


dR 


dC 


da 

da 


^da 


^da^ 


da 


(74) 




da 


^dR^^dr) dR~^ 


da 


^dR^^ dr) dt 


and by putting the coefficients of {a^,R^,T^,aR,dT,TR} and constant terms, 
we obtain the following system of partial differential equations for X components 

{a,/ 3 , 7 }: 


: afn + jdafuR + ya/ijT + ^aUafn + o^fRuPa + laa^fRT = 0, (75) 

IrRO-R = 0, fRTOCT = 0, aafRR = 0 (76) 

dR : u/rrP + a'yfrTT + afRRaa + apRfRR (77) 

+aiRfRT + “^aRfR + 2 /rru = 0 , 

dT : afRRRP + afRRTj + (2 a + aaa + uPr) Jrr + ajRfRT + 2aRfR = 0 (78) 
RT : urJrt + arpRR = 0 (79) 

(jR +P{T-ip{a) + p{a))^afRT (80) 

+70 (t- 3p (a) + p (a)^ Jtt + PaRfRR 

+3(^{T-3p{a) + p{a) + ^ap'(a) - ap'{a)^fT - f + ^ap'(a) + p{a) + R/ru = 0. 


It is a hard job to find all possible solutions of this nonlinear system of first order 
coupled partial differential equations (PDEs). Inspired directly from the case of 
general relativity we’ll limit ourselves to the following simple cases: 

• The case of Einstein gravity with matter components /(i?, T) = i? + 2A + 
g{T): because any theory of modified gravity must be reduced to Einstein- 
Hilbert action at low curvature regime, we are interested to study a solution 
of f{R,T) in which the Einstein-Hilbert term is dominated as the leading 
term of theory. If we substitute f{R,T) = R + 2A + g{T), we obtain the 
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following exact solution: 


p 

II 

II 

-q 

II 

0 

(81) 

g{T) = ^ 

(82) 

p{a) =2K+^. 

(83) 


We mention here that this is an special class of solutions founded before 
in literature [SI] (for the case of purely f{R,T) = F{R) see |ST]). The 
associated Noether charge reads: 


— l‘2Cld^/(i — Sq- 


(84) 


The corresponding scale factor is obtained as follows: 




1 12/3 


1 ^Hot 


Hn = 


SCi 


This solution can be written in terms of q-exponential family 

34 

a{t) = aoe 2 / 3 (Hot), oq = (-^)^^^- 


(85) 


( 86 ) 


The pressure term deserves more investigations. The first term implies on 
the existence of the ’’background pressure” po = 2A, which it can be realized 
by expectation vacuum energy of some quantum fields. The second term 
is a dark matter term, if we identify ci = Pmo, as the dark matter density 
at present era, t = 0, oq = 1. In the case of perfect fluid with equation of 
state w = we have p ^ By a power-law expansion, a{t) ^ F, 

to be accelerated Universe we must have p > 1, in our case the solution is 
not accelerating solution. 

• Case with f{R,T) = h{R) + k{T): if we would like to pass to modified 
gravity, we should try to solve systems of equations by assuming that R 
h{R), and by taking into the account the matter sector k{T). This simple 
assumption gives us the following exact solutions for the model m under 
Noether symmetry approach: 

a = 0, 13 = ^^, 7=-A(r), (87) 

a a 

f dT 

h{R) = Cl I {R - T)Fi(r)dr, k(T) = C, J ^ + C^. (88) 


here {A(T), Fi(i?)} stands for a pair of arbitrary functions. The conserved 
Noether charge associated with this model reads as follows: 


cro = -6CiaaFi{Rf. 


(89) 


Due to the leakage of more information of the form of Fi{x), we cannot 
integrate it. But it is possible to solve it for Starobinsky inflationary model 
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F{R) = R + solution as the following: 


t + C2 = 


fa{t) 


xdx 


lo R{-x^-2Aj{y))+4Ci 


here 


R{-x^ - 24j(y)) = Root(^ - x^ - 24j(y) = 0\y=x'^, 
and we define 


(90) 

(91) 


J\y) — 2J 2 A/TR/^i 


aoM^Ci^ 36 

1 ao\/9Ci'^y‘^M‘^ - Ciy^ _ v^g p 

648 


M^Ci^ n664:CfM<^ 

+V9M-^C^\/9Cfym2 - Ciyao) + 

V rTri 


In 


M^Ci 
r{-^Ciao + 9y^C^M^) 
K 9 CiM 


6. Noether symmetry for mimetic /(i?) gravity 


In this section we’ll see how Noether symmetry gives us useful information about 
mimetic f{R) theory. Especially we would like to search for the possible forms 
of potential function V{(f)) in this model. Thanks to the normalization condition 
dy,(j)d^(j) = 1, the potential is an implicit function of t. So, generally speak in the 
dynamical system Lagrangian is time dependent. Let us to start by the same method 
as we used in the previous section. We write the following condition of dZID for f{R) 
mimetic model: 


dC ^dC dC f .da ^da ; 9a \ dC 

f.dp -dp ■dP\dC f.d-i -97 ■d-f\dC „ 

V da dR d(l)J dR ^ da dR d(l)J d^ 


(93) 


where for simplicity we assume that A = X(t),p = p{a). Using (I33L we can write 
the following system of differential equations , linear in {a,/3,7}: 


—SaRfii — PaRfim + 3a/ — 3aA — 3aU (p) — — aap'{a) — Sap — 0(94) 


3a + 2a70 = 0 (95) 

+ pafRR + oPPafRR + afn = 0 (96) 

ckr/rr = 0 (97) 

Aa^7o - Qa^fR - SaP^fuR = 0 (98) 

2aaRfR + 2aafRR + paP /rrr + a^aafRR + a^PRpRR = 0 (99) 

XajR — Sa^fRR = 0 . ( 100 ) 


which are obtained setting to zero the coefficients of the different terms 
{a^, dR, dp, pR}. We can distinguish some possible cases: 
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• Case of GR: the possible non trivial solution for f{R) = i? + 2A read as 
the following: 


v{4>) = c^ 


a (a, R,(/>) = - Cj sj ciA sin 

a( r> J.\ (2 Cl + 3) Cj a'^i^/ciA / 

P (a, R, (j)) = - — -I 


7 (a, R, (j)) = Cl I C 2 sin 


. f 

sin -:-( 


V 2 


p{a) = 2 A — A — C 5 . 


( 101 ) 

_ C 2 cos > (102) 

C 3 sin - O 2 cos 

+ Cs cos > (104) 

(105) 


This solution corresponds to a constant potential form. The Noether con¬ 
served charge reads as the following: 

Q = Ci \/ciA ^€3 sin ~ cos ^^106) 


+ 2 \Ci a'i+3 ( C 2 sin 


. ( 


2 

-|- C 3 cos 


( 


Because ^ = t the equation is integrable to give us a{t). For Q = 0, the 
exact solution existed for the scale factor: 


a{t) = ail {C 3 sin (wt) 


C 2 cos {ujt)y 


Cl 


V Cl A 

w = . (107) 


It defines an oscillatory solution with Type IV future singularity [53]. If 
C 2 = 0, w —>■ iw, n = |, It is identified in the late-time ACDM era. 
Following [^, we can classify the future singularities as follow: 

— Typel: (’’BigRip”): f —>• ts,a —>• 00, p — >■ 00 and |p| —>■ 00. 

— Typell: ("sudden”): t —>■ G,a as,p ^ Ps and |p| —>■ 00. 

— Type III : t —>■ GjO —>■ Us,p —>■ 00 and |p| —>■ 00 

— Type IV : t —ts,a —>■ Os, p —>■ 0 and \p\ —>■ 0 and higher derivatives of 
H diverge. 

Here tsjOs andps are constants with as ^ 0. 

For our scale factor (ITOTI) the Hubble parameter and first and second deriva¬ 
tives of R read as the following: 


H{t) = — = nuj 
a 


C 3 + C 2 tan(a;t) 


-C 3 tan(a;t) — C 2 - 
H = —nuP{C2 + C^) sec^(wt)(Ca tan(a;t) — C 2 )~'^ 

2nui^{C2^ + C'3^)(sin(wf)C'2 -|- Cscos^uit)) 


H = 


Cf sin^(a;t) — Cf cos^iwt) — 3 ( 73(72 sin(wt) — C 2 cos(wt) 


(108) 

(109) 

110 ) 
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There is solution(s) for H{ts) = oo where 


= nTT ± arctan 


<^ 2,3 

Lc-l + clJ’ 


n G Z. 


Where 


(111) 


lim ^ /JJ2) 


So, our solution represents Type IV future singularities. It is remarkable 
that scale factor remains finite at t —>■ Also, because we must have 

limp(a)= lim ( 2 A — A — Cb + = 0 (113) 

t—^ts a—^as \ / 

we obtain: 


2 A-A-+ 0/(75 = 0 (114) 

the cosmological constant term is ’’calibrated” as the following: 

A= i(A+Q-o/Cs). (115) 

• Modified gravity : if we put f{R) = R + 2A + h{R) a very careful analysis 
of the system of equations gives us the following solutions: 


6.1. Solution with quadratic potential form: Hybrid 
inflation model 

The system of equtions has the following exact solutions for a set of func¬ 
tions: 


V (</>) = C5 + Ce , 

(116) 

2 

a{a,<j)) = --aCi, 

(117) 

l3{a,R,(j)) =2-^, 

(118) 

"f{a,R,(j)) = Ci(j)+ C 2 , 

(119) 

h{R) = CsR+C^, 

(120) 

p (a) = — C 5 + 2 A + — A + Cj. 

(121) 


The quadratic potential is for a massive scalar. This type of potential used 

as an inflationary model, so called as a Hybrid inflation model. More pre- 

2 

cisely, if we set C 2 = 0, = Vo,Ce = the potential in this regime is 

written as: 

V{^) = Vo + 


(122) 
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This case is a vacuum energy, drive inflation, and is distinct from the other 
classes, like large-field or small-field models [65]. Conserved charge is writ¬ 
ten as the following: 

S* = SCiCaa^a -F 2\a^{Cit + C 2 ). (123) 


Specially by integrating this equation for E* = 0 we obtain: 

_ A t(Ci t+2 C 2 ) 

a(t) = age . 

This model is described as the bouncing universe [66167) a(t) ^ 
show this equivalence, we rewrite it in the following form: 

a(t) = dge ^ ao = aoe^V?, 


(124) 
To 

(125) 


Hubble parameter and deceleration parameter read as the following: 


H = - 


xcl 

4C3Cf 


t, '? = -(! + 




(126) 


To have acceleration expansion we should have: 

ACf ACfty „ 
“ ACsCf “ 4^3 ^ 


(127) 


We obtain the following cases: 

— If ^ > 0 then we have acceleration for time intervals t € 

— If ^ < 0 then we have acceleration for time intervals t € 

In the above discussions, we only consider the case that the universe is 
expanding. 


6.2. Solution with exponential inflationary models 
A possible solution is the following hyperbolic model: 

V ((/>) = C6+2Cy cosh 

a{a,(j)) = -T \fC 3 smk , 


7 (a, i?, (/)) = - ia'^^^/Cs cosh^^^-^-^, 
^ (a, R, ^) = sinh 


3CiR 


Cl 


h{R) = -{jXC^ Ci^ + 1)R+C5, 
p (u) = — (A-|- Cg -\- 6A -|- 6 Cy -|- 3 Gg ) -|- CgO^ 


(128) 

(129) 

(130) 

(131) 

(132) 

(133) 
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Because all functions must be real valued, so it is adequate to write C 3 = 
—M^, where M is a free parameter. The associated conserved charge reads: 


Qo = a 


— „2+C4 


.2(4-C'4C'i 2)M . ^ 2(t + C2) ,,,, 2(t + C2\ 

a— -— smh ^ ' + 2MXa cosh ^ 


Cl Cl C: 

It is completely integrable,specially when Qo = 0, we obtain: 


^|34) 
1 -I 


lit) = Cl (^sinh(^ 2 ^^^^ , 


ACi^ 


2(-4 + C4Ci") 


(135) 


In case of n = | it coincides with the late-time ACDM era. 


7. Conclusions 

Motivated by recent observational data, indicates that we live in an accelerating 
universe, several forms of modified gravities have been proposed to resolve and ex¬ 
plain this physical phenomena. One of the most popular and physically acceptable 
candidates is f{R) gravity and its extensions. In our work we established Noether 
symmetry issue for two types of f{R) theories: a type of non-minimally coupled 
model is called as f{R,T) and mimetic /(R). We started by reviewing the ba¬ 
sic physical foundations of these theories. In f{R,T) model we have been written 
point-like Lagrangian for flat FLRW metric. We studied equations of motion and 
Noether symmetry form for it. Two important classes of solution for /(R, T) were 
found. In the first class, we show that the generalized q-exponential scale factor is 
an exact solution which it mimics the background with the background pressure. 
Other solutions were found as general family of additive models, with an exact so¬ 
lution for scale factor in terms of elementary functions. This cosmological solution 
was obtained by considering the Starobinsky model /(R) = R -I- aR^. In mimetic 
/(R) theory, we have been considered Noether symmetries. We observed that there 
are two classes of solutions: the first is equivalent to the GR with dark matter 
and the calibrated cosmological constant. This case mimics a type of cosmological 
solutions with type IV future singularities, where higher derivatives of H diverge. 
Another case is modified gravity with two specified forms of potential functions: 
hybrid inflationary model,in which the scale factor evolves in the bouncing scenario. 
The second family is exponential form, in this case scale factor mimics the form 
of ACDM model perfectly. So, all cosmological models including ACRM,bouncing 
and oscillatory solutions with future singularities are described perfectly by Noether 
symmetrized /(R, T) and mimetic /(R) theories. We conclude that Noether sym¬ 
metry is able to provide a very excellent way to study cosmological implications of 
extended /(R) theories. 
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